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I.  INTRODUCTION 


Optical  wavefronts  are  aberrated  as  they  propagate  through  the  turbulent 
atmosphere.  These  aberrations  limit  the  performance  of  large  aperture  ground 
based  optical  systems.  Adaptive  optical  systems,  which  measure  and  correct 
the  atmospherical  ly  induced  aberrations  in  real  time,  can  significantly 
improve  the  performance  of  ground-based ,  large-a  perture  optical  systems.  The 
key  components  of  an  adaptive  optical  system  are  the  wavefront  sensor,  which 
measures  the  aberrations,  and  the  wavefront  corrector,  (e.g.,  deformable 
mirror)  the  figure  of  which  is  adjusted  to  compensate  for  the  aberrations 
induced  by  the  atmosphere.  Since  the  earliest  days  of  adaptive  optics 
research  there  has  been  debate  over  whether  the  desired  figure  on  the 
wavefront  corrector  should  be  implemented  with  nodal  (i.e.,  local)  or  modal 
( i  .e .,  global  )  influence  functions.  The  objective  of  this  paper  is  to 
compare  nodal  and  modal  wavefront  correction. 

In  the  next  section  a  discussion  of  the  attributes  of  Zernike  polynomials 
is  presented.  In  the  third  section  we  review  the  concepts  of  nodal  versus 
nodal  correction.  Following  that,  the  analyses  used  to  compare  the 
performance  of  nodal  and  modal  wavefront  correctors  are  presented.  The 
results  are  presented  in  section  5  in  the  form  of  plots  which  compare  the 
performance  of  the  two  approaches  as  a  function  of  various  parameters. 

We  conclude  with  observations  on  the  feasibility  of  implementing  either 
nodal  or  modal  correction  for  cases  of  interest. 

II.  ATTRIBUTES  OF  NODAL  AND  MODAL  CORRECTION 

From  a  theoretical  aspect,  nodal  and  modal  correction  are  very  similiar. 
The  wavefront  error,  (fe  for  either  approach  can  be  written  as: 


0c(r)  is  the  implemented  wavefront  correction. 

To  implement  a  nodal  correction,  the  wavefront  to  be  corrected  is 
decomposed  using  a  basis  which  is  determined  by  the  nodal  (actuator) 
influence  function  of  the  wavefront  corrector.  This  decomposition  results 
in  a  set  of  coefficients  which  correspond  to  the  drive  signal  required 
at  the  corresponding  node  to  implement  the  desired  correction.  The 
coefficients  are  typically  chosen  to  minimize  the  mean  square  error  between 
the  aberrated  wavefront  and  the  implemented  correction. 

For  modal  correction,  the  wavefront  is  decomposed  into  various  modes.  In 
this  paper  the  Zernike  polynomials  will  be  used  as  the  basis  for  the 
decomposi tion .  As  for  the  nodal  case,  the  decomposition  results  in  a  set 
or  coefficients  such  that  each  coefficient  corresponds  to  the  strength  of  a 
particular  node  which  is  contained  in  the  aberrated  wavefront.  As  before, 
the  coefficients  are  typically  chosen  to  minimize  the  mean  square  error 
between  the  aberrated  wavefront  and  the  implemented  correction. 


Hence,  for  both  nodal  and  modal  correction,  the  following  equation  for 

2 

mean  square  error  averaged  over  the  aperture,  (T  ,  is  used  to  measure  how 


good  the  correction  is: 
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where:  W(r)  is  a  window  function  which  defines  the  aperture  of 
the  wavefront  corrector,  and 


~> 

W(r )  = 


1  inside  the  aperture 
0  outside  the  aperture 


r 
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the  <  >  brackets  indicate  an  ensemble  average. 

A  mathematical  description  of  a  nodal  correction  is  given  below: 


viv-cv-; 


>  -> 


t(r)  =  a  p (  r  -  r  ) 


where:  a  is  the  drive  signal  applied  at  location  r 

i  i 

and,  p(r)  is  the  influence  function  of  the 

wavefront  corrector.  (Implicit  here  is  the 

assumption  that  the  influence  function  is  the 

same  for  every  node,  which  is  generally  the 

case  ) . 

Usinq  equation  (4)  in  equation  (1)  and  the  result  in  equation  (2)  gives 

2 

the  residual  mean  square  error  for  nodal  correction,  d  . 

Ne 


r  ~>  •*>  2  -> 

</w(  r  )  [  (J(r)  -  a  p(  r  -  r  )  ]  dr  > 

i  =1 


The  error  given  by  equation  (5)  can  be  set  to  any  desired  level  by 
selecting  the  proper  number  of  nodes.  In  general,  the  greater  the  number 
of  nodes,  the  lower  the  residual  error  will  be  for  a  particular  W(r). 
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The  relationship  between  residual  error  and  the  number  of  nodes  will  be  put 
on  a  sound  mathematical  framework  in  section  4. 


The  residual  mean  square  error  for  a  modal 
below: 


2 

corrector,  (T  ,  is 
Me 


shown 


2 
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Me 
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->  2  -> 

b  Z  (  r  )]  dr  > 
j  J 
j=l 


</w(  r  )[  0  (r 


(6) 


where:  Z  (r)  is  the  jth  mode  Zernike  polynomial  (defined  below) 

j 

and  b  is  the  strength  of  the  jth  mode. 

j 


In  general,  the  greater  the  number  of  modes  which  are  corrected,  the  smaller 
the  residual  error  will  be.  The  desired  residual  error  can  be  achieved  by  proper 
selection  of  the  modes  which  are  implemented  in  the  wavefront  corrector.  This 
relationship  will  be  presented  in  a  mathematical  form  in  section  4. 


As  shown  in  equations  5  and  6,  the  theoretical  performance  of  both  nodal 
and  modal  correctors  is  the  same  (i.e.,  any  level  of  residual  error  can 
be  obtained  if  enough  nodes  or  modes  are  used  in  the  wavefront  corrector). 

The  essence  of  the  debate  between  the  use  of  nodal  or  modal  correction  is 
thus  the  ease  with  which  the  desired  correction  can  be  implemented. 


Conceptual  one  dimensional  nodal  and  modal  wavefront  correctors  are  shown 


in  figures  2-1  and  2-2.  Nodal  correction  implies  a  local  deformation  of  the 
surface  of  the  corrector  as  a  result  of  the  application  of  a  drive  signal; 
modal  correction  implies  a  global  deformation  of  the  surface  as  a  result  of 


the  application  of  a  drive  signal.  The  conjecture  which  is  often  made  is 
that  the  number  of  drive  signals  required  for  a  modal  corrector  could  be 
less  than  the  number  of  drive  signals  required  for  a  nodal  corrector.  In 
the  following,  we  investigate  whether,  for  atmospherically  induced  aberrations 
it  is  likely  that  modal  correction  will  require  fewer  drive  signals  than  a 
nodal  corrector. 


figure  2-1.  Conceptual  Nodal  Wavefront  Corrector 


Figure  2-2.  Conceptual  Modal  Wavefront  Corrector 


III.  PROPERTIES  OF  ZERNIKE  POLYNOMIALS 


Zernike  polynomia  1  s ,  or  Zernike  radial  polynomials,  are  a  set  of 
polynomials  which  are  defined  on  a  unit  circle.  They  are  a  special  case  of 
•Jacobi  hype  rgeome  tri  c  polynomials  which  are  used  to  establish  a  base  set 
of  two~d i men siona  1  polynomials  that  can  form  a  polynomial  function  which 
is  a  product  of  both  a  purely  radial  and  purely  angular  function  [1], 


The  purely  radial  function,  R(r),  describes  how  the  surface  changes  from 


the  center  to  the  edge  of  the  circular  aperture  and  the  purely  angular 
function,  Q{»),  describes  how  the  surface  changes  while  moving  azlr.jthally 
around  the  aperture.  Since  most  optical  systems  have  circular  apertures, 
Zernike  pol ^lomial  s  provide  an  effective  way  of  expressing  phase 
aberration  over  these  apertures. 


Zernike  polynomials  are  conveniently  expressed  in  polar  coordinates  and 
m 

can  be  symbolized  as  R  (r)  where  n  is  the  highest  power  of  r,  or  radial 

n 

degree  and  m  is  the  azimuthal  frequency,  or  angular  order  for  which  the 
angular  function  repeats  itself.  Table  3-1  shows  the  first  19  Zernike 
polynomials  expressed  in  both  polar  and  cartesian  coordinate  systems  [2], 


The  polynomials  are  defined  here  by  [3]: 

m 


a)  Zeven  =  4n  +  1  R  (r)j2  cos  (mo) 
j  n 


- i  n 

b)  Zodd  =  Jn  +  1  R  ( r ) 42  sin(mo) 


Z  =  Jn  +  1  R  (r ) 
j  n 


-i  o 


m  =  0 


m  =  0 


(1 


(n-m)/2  s 

m  (-1)  (n  "  s)! 

whe  re :  R  ( r )  =  y  -  - 

n  ^  s!  [(n-Hn)/2  -  s ] !  [(n-m)/2  -  s]! 

s=0 


n-2s 

r  (! 


The  index  j  is  a  mode  ordering  number  and  the  values  of  n  and  m  are  always 
integers  which  satisfy  m  <  n,  and  n  ~m  =even.  Consequently,  only 
polynomials  with  certain  combinations  of  n  and  m  exist  (see  table  3-2)  [3], 


It2-! 

2x2+2y2-l 

r2  cos20 

7^-J2 

r2  sin20 

2xy 

(3r2-2)r  cos0 

3x3+3xy2-2x 

(3r2-2)r  sin0 

3x2y+3y3-2y 

6r4-6r2+l 

6x4+12x2y2  +6y4-6x2-6y2+l 

r3  cos30 

x3-3xy2 

r3  sin3t 

3x2y-y3 

(4r2-3)r2  cos20 

4x4-3x2+3y2-4y4 

(4r2-3)r2  sin20 

8x3y+8xy3-6xy 

(10r4-12r2+3)r  cos0 

10xs+20x3y2+10xy4-12x3-l2xy2+3x 

(10r4-12r2+3)r  sin0 

10x4y+20x2y3+10y5-12x2y-12y3+3y 

20r6-30r4+12r2-l 

20x6+60x4y2+60x2y4+20y6-30x4 

-60x2y2-30y4+12x2+12y2-l 

• 

r4  cos40 

x4-6x2y2+y4 

r4  sin40 

4x3y-4xy3 

(5r2-4)r3  cos3 

5xs-10x3y2-4x3-15xy4+12xy2 

(5r2-4)r3  sin3 

-5y^-10x2y3+4y3+15x4y-12x2y 

Table  3-1 


As  mentioned  earlier,  Zernike  polynomials  are  particularly  valuable  for 
their  unique  properties  over  a  circular  aperture.  Also,  they  are  extremely 
valuable  due  to  their  relationship  to  classical  aberrations.  Therefore, 
they  provide  a  convenient  mathematical  expression  of  the  aberrating  content 
in  a  wavefront  using  familiar  terms  (see  figures  3-1  through  3-6)  [4]. 

The  properties  posessed  by  the  Zernike  radial  polynomials  over  a  circular 
aperture  are  enumerated  below  [2]: 

(1)  They  are  orthogonal  over  the  unit  circle. 

(2)  They  are  normalized. 

m 

(3)  In  the  term  R  (r),  if  m  is  even  n  must  be  even, 

n 

and  if  m  is  odd  n  must  be  odd. 
m 

(4)  The  term  R  (r)  contains  no  power  of  r  less  than  m 

n 

(i  .e.  n  >  m). 

As  described  above,  Zernike  polynomials  possess  several  properties  which 
make  them  ideally  suited  for  describing  the  aberrations  present  in  an  optical 
system.  From  (1),  one  of  their  most  important  properties  is  that  they  form 
an  orthogonal  set  of  polynomials  over  a  normalized  circular  aperture.  This 
implies  that  each  term  is  independent  from  all  others.  Consequently,  varying 
one  term  in  the  series  will  not  effect  the  values  of  the  other  terms.  This 
characteristic  is  extremely  important  in  the  process  of  correcting  an 
atmospherically  induced  aberrated  wavefront  in  an  adaptive  optical  system. 
Lastly,  each  term  of  the  set  of  Zernike  polynomials  can  be  related  to  another 
set  of  deformations  of  the  lowest  order  called  the  Seidal  aberrations  [21. 
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Figure  3-4b 


ZERNIKE  POLYNOMIALS 
FIGURES  WITH  COEFFICIENTS 


F  i  gure 

Name 

Figure 

Name 

3-1) 

a . 

Tilt 

3-4)  a. 

Coma 

b. 

b. 

3-2) 

De  focus 

3-5) 

Spheri  cal 

3-3) 

a . 

Asti  gmati sm 

3-6) 

151st  Zernike 

b. 

Polynomial 

Table  3-3 


IV.  MODAL  TO  NODAL  COMPARISON 

In  this  section  we  compare  the  number  of  drive  siginals  required  for  the 
nodal  and  modal  methods  to  implement  a  wavefront  correction  which  equally 
compensates  for  atmospherically  induced  wavefront  deformations. 


The  equation  for  the  mean  square  residual  error  of  the  wavefront  using 
nodal  correction  is  defined  here  by  [5]: 

2  5/3  2 

O'  =  a  (d/rQ  )  [red  ]  (It 

Ne 

2 

where:  0  is  the  mean  square  residual  error  of  the  wavefront. 

Ne 

a  is  a  constant  which  depends  on  the  wavefront  corrector 
(.141  <  a  <  .34) 


d  is  the  spacing  between  actuators  in  the  wavefront  corrector, 
and 

r0  is  the  phase  coherence  length  of  the  atmosphere. 

The  equation  for  the  residual  error  of  the  wavefront  using  modal 
correction  is  defined  here  by  [3]: 

2  -,[3/2  5/3  2 

(T  =  .2 944 J  ( D/r0  )  [rad  ]  ( 11) 

Me 

2 

where:  (T  is  the  mean  square  residual  error  of  the  wavefront. 

Me 

J  is  the  number  of  modes. 

D  is  the  diameter  of  the  aperture. 

ro  is  the  phase  coherence  length  of  the  atmosphere. 

The  two  methods  of  wavefront  correction  were  compared  to  one  another  by 
plotting  various  parameters  of  the  two  equations.  The  series  of  four  graphs 
developed  for  the  comparison  will  be  discussed  here.  The  analysis  and  results 
of  this  graphical  comparison  will  be  presented  in  the  next  section. 

OBSERVATIONS  FROM  GRAPH  4-1 

The  first  graph  (graph  4-1)  is  a  plot  of  the  mean  square  residual  error 

2 

for  modal  correction  versus  t.he  number  of  modes  ( J)  .  (T  was  calculated 

Me 

from  equation  (11)  for  values  of  D/rc  shown. 

This  oraph  shows  the  number  of  modes  (J)  necessary  to  achieve  a  given 
mean  square  residual  error.  For  each  of  the  four  cases  of  D/r  =  10,  20  ,  50  , 

O 

and  100,  the  plots  are  exponential  and  show  that  as  the  number  of  modes  (J) 


increases,  the  residual  error  squared  ( (T  )  decreases  and  conversely,  J 
2  Me 

decreases  as  (T  increases.  Therefore,  the  greater  the  number  of  modes 

Me 

corrected,  the  less  the  wavefront  error  will  be. 


Further  observation  shows  that  as  the  ratio  of  the  diameter  of  the 
aperture  to  the  phase  coherence  length  of  the  atmosphere  (0/ro  )  increases, 
the  number  of  modes  (J)  necessary  to  obtain  an  acceptably  small  wavefront 
error  rapidly  increases.  Consequently,  in  order  to  achieve  the  image  quality 
desired  in  an  optical  system,  the  number  of  mode  corrections  necessary  may 
be  very  high  for  the  D/r  of  useful  systems. 


OBSERVATIONS  FROM  GRAPH  4-2 


The  second  graph  (graph  4-2)  is  a  plot  of  the  mean  square  residual  error 

for  nodal  correction  versus  the  actuator  spacinq  divided  by  the  phase 

2 

coherence  length  of  the  atmosphere  (d/ra  ).  (7  was  calculated  from  equation 

Ne 

10  for  a  =  0.2, 0.4,  and  1.0  at  values  of  d/rc  =  0.2,  0.5,  1.0,  and  2.0. 

Each  of  the  3  plots  contains  4  data  points. 


Each  of  the  three  plots  is  exponential  in  shape  and  indicates  that  as 

d/rQ  increases,  the  residual  error  of  the  wavefront  also  increases.  In 

addition,  analysis  of  the  three  plots  indicates  that  as  the  parameter  (a) 

2 

increases,  the  larger  (7  becomes  for  a  given  value  of  d/ro  .  This  observation 

Ne 

follows  from  the  nodal  equation  (1)  seeing  that  the  term  (a)  acts  as  a 
multiplication  factor. 


Further  observation  of  the  3  plots  suggests  that  the  implementation 
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of  a  nodal  wavefront  corrector  in  a  compensated  imaging  system  can  achieve 

2 

an  acceptable  mean  square  wavefront  error  of  0.6  rad  ,  for  a  relatively  large 
d/i^  (equal  to  approximately  2.0)  and  a  value  of  the  parameter  (a)  equal  to  0.2. 

OBSERVATIONS  FROM  GRAPH  4-3 


The  third  graph  (graph  4-3)  is  a  plot  of  (d/r  )  verses  the  number  of 

2  2  2 

modes  (J)  for  the  case  where  d  =  O’  .  O’  was  calculated  from  equation 

Me  Ne  Ne 


2 

10,  for  a  =  0.2  and  d/r  =  0.2  to  2.0  in  0.2  increments.  O’  was  then 

Ne 

substituted  into  equation  11  with  D/r0  =  10,  20,  50,  and  100  to  calculate 
the  corresponding  value  of  J.  Consequently,  each  of  the  4  plots  contains 
10  data  points.  The  d/re  scale  is  linear  from  0.2  to  2.2  and  the  J  scale 
is  logarithmic  from  1  to  1,000,000.  Graph  3  provides  a  comparison  betwee 
the  performance  of  modal -to -nodal  wavefront  correctors. 


On  the  x  axis  is  the  number  of  modes  (J)  which  was  derived  from  the 

modal  equation  and  on  the  y  axis  is  (d/r0  )  which  was  derived  from  the  nodal 

equation  (10).  Graph  4-3  contains  4  plots,  for  D/re  =  10,  20,  50,  and  100, 

which  are  exponential  in  shape.  From  the  four  curves,  it  is  observed  that 

the  greater  the  number  of  modes  (J),  the  smaller  i/r0  must  be  to  provide 

the  same  performance  of  the  wavefront  corrector.  For  a  residual  error  of 

the  wavefront  specified  by  a  given  value  of  d/r  ,  the  greater  D/r  becomes, 

2 

the  greater  the  number  of  modes  necessary  to  achieve  this  (7  . 

e 


Usinq  the  relationship 

2 

N  =  (  D/d  )  (12) 


where:  N  is  equal  to  the  number  of  nodes. 
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A  direct  comparison  can  be  made  between  the  number  of  nodes  (N)  to  the  number 
of  modes  (J)  necessary  to  achieve  a  given  residual  error  of  the  wavefront. 
Solving  (12)  for  d  gives: 


d  =  0/  J7T 

Using  (13)  in  (10)  gives  : 


-5/  6  5/3 

=  a  N  (D/r0  ) 


Equating  (15)  and  (11)  gives: 


-5/6  5/3  -J3/2  5/3 

a  N  (D/r0  )  =  .2944  J  (D/r0  ) 


solving  for  N  and  simplifing  gives: 


6/5  1.04 

N  =  4.338  a  J 

A  plot  of  this  function  is  shown  in  graph  4-4. 


(13) 


(14) 

(15) 


(17) 


Using  the  above  we  can  finally  compare  the  feasibility  of  using  nodal 

verses  modal  correction  for  atmospheric  turbulence.  To  make  the  comparison  we 

use  a  hypothetical  adaptive  optical  system  which  has  a  diameter,  D,  of  2 

meters.  We  assume,  based  on  measurements,  an  atmospheric  coherence  length, 

of  10  cm.  Based  on  measured  data,  a  reasonable  value  for  the  parameter  (a  )  i  s 
2  2 

0.2.  We  select  O’  =  0.561  rad  as  a  reasonable  residual  wavefront  error. 

Ne 


2 


Using  the  above  values  for  D,  r  ,(T  and  a  in  (15)  gives  an  N  of  116. 


d/ri  VALUES 


The  spacing  between  the  actuators  in  the  wavefront  corrector 
divided  by  the  phase  coherence  length  of  the  atmosphere  as  a 
function  of  the  number  of  modes 


J  VALUES 

Graph  4-3 


The  number  of  nodes  as  a  function  of  the  number  of  modes 
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Substituting  this  value  of  N  into  (17)  gives  151  for  a  value  of  J.  A  plot 
of  this  Zernike  polynomial  is  shown  in  figure  3-6.  It  would  appear  that 
many  actuators  would  be  required  to  implement  this  mode.  The  data  from  the 
comparison  done  on  graph  4-3  is  shown  in  table  4-1. 

NUMBER  OF  NODES  TO  NUMBER  OF  MODES  COMPARISON 


10 

20 

50 

100 

2,239 

10,036 

63,673 

258,453 

2,600 

11,000 

65,000 

250,000 

107 

386 

2,239 

10,036 

140 

480 

2,600 

11,000 

25 

107 

606 

2,239 

35 

140 

740 

2,600 

Table  4-1 

From  table  4-1  it  can  be  seen  that  the  number  of  nodes  (N)  increases  for  a 
given  d/rQ  and  larger  values  of  D/  rQ  .  This  was  also  the  case  for  the  number  of 
modes  J  as  was  stated  earlier.  Overall,  the  number  of  modes  is  approximately 
the  same  as  the  number  of  nodes  for  producing  a  particular  residual  error 
of  the  wavefront  for  all  combinations  of  d/r^  and  0/ rQ  .  Consequently,  an 
atmospheric  modal  wavefront  corrector  will  require  more  actuators  than  a  nodal 
corrector  to  compensate  for  a  heavily  aberrated  wavefront  in  a  large  aperture 
ground-based  optical  system,  since  each  mode  will  require  several  actuators. 


V.  CONCLUSION 

In  theory  modal  correction  systems  are  simpler  to  implement  than  nodal 
correction  systems.  For  very  low  mode  aberrations,  (e  .g.  ,ti  It  .focus )  modal 
corrections  show  that  this  simplification  can  be  realized.  However,  for 
large  aperature,  visible  wavelength,  adaptive  optical  imaging  systems  the 
number  of  modes  which  are  required  to  be  corrected  becomes  very  high. 
Implementation  of  each  of  these  modes  with  only  a  few  actuators  does  not  seem 
practical.  Therefore,  for  such  systems,  nodal  wavefront  correctors  are  the 
wavefront  corrector  of  choice. 
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